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Abstract. Crosslinks and molecular motors play an important role in the 
organization of cytoskeletal filament networks. Here we incorporate the effect of 
crosslinks into our model of polar motor-filament organization [Phys. Rev. E 71, 
050901 (2005)], through suppressing the relative sliding of filaments in the course 
of motor- mediated alignment. We show that this modification leads to a nontrivial 
macroscopic behavior, namely the oriented state exhibits a transverse instability in 
contrast to the isotropic instability that occurs without crosslinks. This transverse 
instability leads to the formation of dense extended bundles of oriented filaments, 
similar to recently observed structures in actomyosin. This model also can be applied 
to situations with two oppositely directed motor species or motors with different 
processing speeds. 
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1. Introduction 

Biological cells consist to a large degree of a complex, self-organizing viscoelastic fluid, 
the cytosol. Its main constituents include cytoskeletal proteins such as actin and tubulin, 
which exist mainly in the polymerized form as semi-flexible actin filaments and stiff 
microtubules [U [2] . The entangled networks of microtubules and actin filaments form 
the cytoskeleton of most eukaryotic cells, stabilize their morphology and determine 
the rheological properties of the cytosol. These intricate networks are created and 
maintained by an efficient mechanism which involves various types of motor proteins 
as well as passive crosslinking proteins. Motors are specialized protein molecules that 
move along the cytoskeletal polymer scaffold and perform directed intracellular transport 
[3]. Additionally, if motors bind to more than one filament, they are able to move 
the filaments and reorganize the cytoskeleton itself. The crosslinks connect different 
filaments but do not move along them. Their main function is thus believed to provide 
rigidity and elasticity to the cytoskeleton. 

Various experiments have been performed in recent years that shed light on the 
viscoelastic behavior of entangled cytoskeletal filament solutions, ranging from filament- 
motor mixtures [1] , crosslinked filaments [HI [6[ [7] , and systems of filaments, motors and 
crosslinks [8]. Surprising new effects have been found such as active fiuidization of actin 
gels by myosin motors [9]. 

Maintained in a state far from equilibrium, the active filaments exhibit a strong 
tendency towards self-organization. Bundles and contracting states have been found in 
vitro in actomyosin extracted from muscle cells [10], and various patterns like ray-like 
asters, spindle-like structures and rotating vortices have been reported in quasi two- 
dimensional mixtures of microtubules and motors [HI [12]. These dissipative structures 
have inspired many theoretical efforts [TSj [HI [151 [HI [13 [HI [H] directed towards 
modeling active filament solutions. 

While crosslinks so far have been mainly investigated only in the context of rheology, 
recently their influence on the dynamics and self-organization also attracted attention 
[20] . In particular, it was shown that crosslinks facilitate the formation of bundles 
in the actin-myosin system: at high concentration of adenosine triphosphate (ATP), 
actin-myosin systems display an isotropic phase; in the course of depletion of ATP 
however, myosin motors become static crosslinks and initiate the formation of oriented 
bundles and cluster-like patterns. Reintroduction of ATP in the bundled state resulted 
in consequent dissolution of the structures and reestablishment of the isotropic state. 

Motivated by these experimental results, we focus here on the effects of static 
crosslinks on the self-organization of polar filaments and generalize the model for 
microtubule-motor interaction introduced in Refs. [T71 [19]. In that model, the 
complicated process of filament interaction via multi-headed molecular motors was 
approximated by instant binary "inelastic collisions", leading to alignment of the 
filament orientation vectors and attraction between their centers of mass. Crosslinks 
alter these interaction rules. In particular, if two parallel filaments are cross-linked. 
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they are not able to slide past each other and become collocated. We model this effect 
here by suppressing relative sliding of the filaments in the course of alignment. Our 
analysis shows that this relatively minor modification produces a nontrivial macroscopic 
effect, namely the isotropic density instability of the polar oriented state of the filaments 
becomes transverse. In the nonlinear regime, this new kind of instability leads to the 
formation of dense oriented bundles, similar to those seen in experiments [20]. In 
contrast, the model without crosslinks demonstrates an isotropic instability in which 
density and orientation of the filaments are uncorrelated, and no bundling occurs. 



2. Model 



Here we outline the model of self-organization of microtubule-motor mixtures developed 
in our earlier works, Refs. [19]. The microtubules are modeled as identical rigid 
polar rods of length L, and the molecular motors are introduced implicitly through 
corresponding interaction probabilities. Binary interactions of microtubules via multi- 
headed molecular motors are approximated by instant inelastic collisions leading to 
alignment of the microtubule orientation angles 0i 2 (or, equivalently, the unit vectors 
^1,2 = (cos 01,2, sin 01^2)) according to the following rules: 

(1) 

Here 0^2 orientations of the two rods after the collision, and the constant 

"restitution" parameter 7 characterizes the inelasticity of the collision (in analogy to 
the restitution coefficient in granular media). The angle between the two rods is reduced 
after the collision by the "inelasticity" factor e = 27 — 1. Of special interest is the totally 
inelastic collision corresponding to 7 = 1/2 or e = 0. In this case the rods acquire the 
same orientation along the bisector n = (cos 0, sin 0), and their center of mass positions, 
ri 2, also align: 











[ ^2 j 










2 



(2) 



r? = r^ = r- = ^:i±^. (3) 

Here 0° 2 and r J 2 are the orientation angles and the center of mass positions after the 
collision. We assume that the alignment through inelastic interaction occurs only if the 
initial angle difference |0i — 02| is smaller than a certain maximum interaction angle 0o. 
For 101 — 02 1 > 00, the angles and the positions are unchanged. The analysis of Refs. 
[T71[T9] showed that in the spatially homogeneous case, the rods exhibited a spontaneous 
orientation transition if the density of the motors (or of the filaments) exceeded a critical 
density. Furthermore, for even higher densities, another instability was predicted which 
is isotropic and leads to inhomogeneous density variations. 

The dynamics of this model can be described by the master equation for the 
probability distribution function P(r,0, t) to find a rod at position r with orientation 
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dt 902 

The first two terms on the right hand side describe rotational and translational diffusion, 
with an anisotropic diffusion matrix of the form 

1 



A, 



[DiiniUj + Dj_{6ij - UiUj)] . 



(5) 



The rotational, Dr, parallel, D\\, and perpendicular, D^, diffusion coefficients for rigid 
rods in a viscous fluid are well known [21]. The third term in Eq. (jll) is the collision 
integral, 

r r r<f>o 

I{r,(f),t)= / dridr2 d4)id4)2W{ri - r2,ni,n2)P{ri, (f)i)P{r2, 

J J J —^Q 

X [5(0 - 0^)5 (r - r1) - 5(0 - 0i)5 (r - n)] , (6) 



where the localization of spatial interactions is introduced through a certain probabilistic 
kernel W{ri - ra, ni, na) [TTl fT9] . 

The kernel W, expressing the probability of interaction between the rods as a 
function of the distance between their midpoints and their orientations, can be obtained 
from the following conditions: (i) since the size of motors is small compared to the 
length of filaments, two rods interact only if they intersect; (ii) due to translational and 
rotational invariance, the kernel depends only on differences 0i —02 and ri — (iii) the 
kernel is invariant with respect to permutations rii n2, Vi —* r2. The kernel can be 
represented as a product of two parts: a part Wq which accounts for spatial localization 
due to the overlap condition of the filaments, and a part describing the motor-induced 
collision anisotropy. 

The first part can be derived from the intersection condition between two rods with 
orientations ni^2- It is easy to verify that the rods overlap if 

|(ri - ra) X nil < L|ni X n2|/2, (7) 

\{ri - r2) X n2\ < L\ni X n2\/2 (8) 

holds. This overlap condition can be expressed in terms of discontinuous G-functions, 

Wo = Wne{L\ni X nsl - 2\{r, - ra) x ni|)e(L|ni x - 2|(ri - ra) x nal) , (9) 

where Wn is a normalization constant, so that j W^dr = 1. Since this discontinuous 
kernel is difficult for calculations, the 0-functions can be approximated by smooth 
Gaussians yielding 



Wo{ri - rs, ni, na) ~ exp 
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where 6 is a cutoff length of order b ^ L. It is convenient to transform the kernel to the 
following representation (the integral of the kernel is normalized to 1): 



Wo{ri -r2,ij) 



Tih^ sin ijj 



exp 
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W^o(R,^), (11) 
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where ■?/' = 0i — 02 is the difference of the orientation angles, and Ry = (ri — ■ n 
and R_L = — (ri — x n are two vectors parallel and perpendicular to the bisector 
direction n. The cutoff length b introduced above can be estimated, for example, by 
comparison of the characteristic kernel width J H^WolH^dTi for the kernels given by 
Eqs. ([9]) and ( fTOl) for some typical angle, say ip = tt/2. Equating both integrals, one 
finds that b^L^ = 2/3. Q 

Finally, the complete kernel can be represented in the form 

W{ri - ra, ni, na) = gWoivi - ra, ^) ( 1 + y (ri - ra) ■ K - na) ) . (12) 



L 

Here g is the interaction rate proportional to the motor density (which can be scaled 
away) and the last term oc j3 describes the anisotropic contribution to the kernel, which 
is associated to the increase of motor density towards the polar end of the filament 
due to dwelling of the motors. Accordingly, the constant f3 can be related to the dwell 
time [IS]. 

Near the threshold of the orientation instability mentioned above, p > pc, the 
master equation can be systematically reduced to equations for the coarse-grained local 
density of filaments p and the coarse-grained local orientation r 

p = J Pir, 0, t)d(P , r = (n) = i- y n P(r, 0, t)# , (13) 
by means of a bifurcation analysis. 



3. Effects of crosslinks in the model 



The effect of crosslinks on the motor-induced interaction of filaments is twofold. First, 
the simultaneous action of a static crosslink, serving as a hinge, and a motor moving 
along both filaments results in a fast and complete alignment of the filaments, as 
shown in Fig. [TJ This justifies the assumption of fully inelastic collisions for the rods' 
interaction. Note that without crosslinks the overall change in the relative orientation of 
the filaments is much smaller: the angle between filaments decreases only by 25-30 % in 
average, see the discussion in Ref. [19]. Complete alignment also can occur for the case 
of simultaneous action of two motors moving in opposite direction, as in experiments 
on kinesin-NCD mixtures reported in Ref. [12], and even for two motors of the same 
type moving in the same direction but with a different speed due to variability of the 
properties and the stochastic character of the motion. Second, the crosslinks inhibit 
relative sliding of rods in the course of alignment, restricting the motion to rotation 
only. Thus, in contrast to the situation considered in Refs. [19] and described by 
Eq. ([3]), in the presence of a crosslink the midpoints of the rods will not coincide after 
the interaction. In fact, the distances 5*1^2 from the midpoints to the crosslink point do 
not change, as it is shown in Fig. [H 

I In our previous works |171 I19j we used a somewhat simpler expression for the kernel, Wq 
exp[— |ri — r2p/6^]. As we have verified, this simplified approximation did not change the results 
on a qualitative level, affecting only numerical prefactors of some nonlinear terms. 




Figure 1. Sketch of the interaction between two filaments, a crossUnk and a molecular 
motor. After the interaction, the motor (shown as a red sphere moving in the direction 
of the red arrow) aligns the filaments along the bisector fi, but the midpoint positions 
do not coincide due to the crosslink (blue sphere). 



To describe the interaction rules in the presence of a crosshnk, we express the 
radius- vector of an arbitrary point on a filament Rj via the position of its midpoint r^, 
the filament orientation rij, and the distance from the center of mass S, Rj = riiSi + rj. 
The intersection point of two rods is given by the condition R* = Ri = R2, which fixes 
the values of 5*1 _2 to 

g,2= ^"^"^^^^^^^ (14) 

Til X '^'2 

Due to the crosslink, the values of do not change during the interaction. Since 
the filaments become oriented along the bisector direction n, the distance of the two 
filament midpoints from the total center of mass will be AS* = {Si — 5*2). Therefore, 
instead of Eqs. ([2D,(I3]) we obtain the interaction rules 

,a .a 7 01 + 02 
01 = 02 = = ^ , (15) 

a ri^r2 . nAS ri + 7-2 . n((ri - r2) • n) 

r'^.. = ±7] = ±1]—^ -. (16) 

2 '2 2 ' 2cosV' 

Here we have introduced the parameter 77 interpolating between two cases: the case with 

crosslinks present corresponds to 77 = 1; for 77 = the previous model, Eqs. ([2]), ([3]), is 

recovered. Thus the value of r] can be roughly interpreted as the effective strength of 

crosslinks or an effective fraction of crosslinks with respect to motors. 

The interaction rules, Eqs. ( fTSjl and ( fT6l) . can be used to evaluate the collision 

integral, Eq.([6]). Omitting lengthy calculations (see the Appendix for details) after 

expanding the master equation (jlj) near the threshold of the orientation instability, we 
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arrive at the following set of nonlinear equations for the coarse-grained density p and 
orientation r: 

= W - CvV - ^ (1 - V V 

- [iSlM + (1 - 2S[MW)^'r' + 2 {S[M - v') ^^^,inT,)] , (17) 

dtT = e{p - Pc)t - A\t\^t + Dr.V^T + Dr,VV ■ T + — ^VV 

47r 



with S[x\ = sin(x)/x and B = b/L. The constants A, e and the critical density pc are 
functions of the maximum interaction angle 0o cind the inelasticity coefficient 7: 

,^2^ (g[0o(27 - 1)] - S[<Po]) {S[Ml + 1)] + S[Ml - 2)] - g(20o) - g(0o)) 
2/00 - (5[20o(7 - 1)] + 5[20o7] - 5[20o] - l)p/27r 

e = ^ [5[0o(7 - 1)] + 5[0o7] - 5(0o) - 1] , Pc = - . (19) 
TT e 

In the following we consider the case 0o = vr as motivated below. Then the density 

equation (fT7|) becomes somewhat simpler since S'[0o] = 0. We have introduced 

rescaled diffusion coefficients, namely Dp = {D\\ + D±)/(2DrL'^) = 1/32, Dt-^ = 

(Dy + 3D^)/{ADrL^) = 5/192 and D^^ = (L)|| - D^)/{2DrL^) = 1/96. In order to scale 

out the motor density g we rescaled density and orientation vector, gp — > p, gr ^ r. 

Also length is normalized by r — > r/L and time by t ^ t/DrL"^. The anisotropic 

contribution proportional to H = pb'^/L'^ = (33"^ is due to the polar distribution of 

the motors along the interacting filaments, while the anisotropic contribution in the 

p-equation is due to the crosslinks. The isotropic higher order diffusion term C^^p 

was included for regularization purposes of the equation at very short wavelengths. 

Assuming additionally 7 = 1/2 (i.e. totally inelastic collisions, as justified above), 

one obtains from Eqs. (fT9|) : e = 4/7r — 1 0.273, A ^ 2.18 and the critical density 

Pc ^ 3.663. 

A sketch of the phase diagram for Eqs. f|T7|) . f|T8l) in the plane of the motor- induced 
anisotropy parameter H and the mean density po is shown in Fig. [21 A uniform isotropic 
state, p = Po and r = 0, loses its stability if po > Pc, independent of the value of H. In 
the spatially uniform case, orientation modulations grow into a polar state with non- 
zero |t| = [e(po — Pc)/AY^'^ and arbitrary orientation of r. Recall that the density 
Po is scaled by the "collision rate" g, and thus is proportional to both the density of 
tubules and the density of motors. This implies that either increasing the number of 
motors or the number of filaments can induce the polar phase. However, in extended 
systems, the growth of spatially inhomogeneous modes leads to the formation of a 
complex state characterized by disordered arrays of vortices or asters depending on the 
value of the anisotropy parameter H [T71 [IS]. Vortices are stable only for small values of 
the anisotropy parameter H; the stability limit of vortices, indicated by the black solid 
line, terminates at a critical point at if = He- The vortex-aster-competition is governed 
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Figure 2. Sketch of the phase diagram of the rescaled density po (being the product 
of motor and filament density) against anisotropy parameter H in the absence of 
crossUnks. Above pc, the polar state is formed. Beyond pd, given by Eq. ([22]). an 
isotropic density instability occurs. Depending on parameters, the density instability 
may happen also prior to the orientation instability, i.e. pd < Pc- In between pc and 
Pd, asters are stable above the critical line (solid black line) while vortices are linearly 
stable below this line. The critical line terminates at the point H — He- In the presence 
of crosslinks, for 0o = bundles occur throughout the polar phases, i.e. beyond the 
red line. However, they are in nonlinear coexistence with the asters/ vortices. 

predominantly by the r-equation, Eq. (fT8|) . and thus prevails whether crosslinks are 
present or not. In the case without crosslinks, ?7 = 0, for densities p > pd, the 
homogeneous oriented state loses its stability with respect to density fluctuations as 
implied by the green dashed line in Fig. [2l If crosslinks are present however, rj = 1, the 
density instability is (to leading order) independent of the value of filament density and 
thus bundles can be found throughout the polar phase, i.e. for all p > pc, where they 
are in complicated nonlinear competition with the aster and vortex defects. 

4. Instability of the homogeneous polar state 

For r] = 0, Eqs. f|T7j) . (fTSj) reduce to the model without crosslinks studied in Refs. p!7l[T9] . 
As it was shown in p!7| [T9] . this equation exhibits an isotropic density instability if 
p > Pd SuS calculated below. In the presence of crosslinks (?7 = 1), the term in Eq. ( ITTIl 
proportional to V^p^ which is responsible for the density instability vanishes, and instead 
a new anisotropic term didj{TiTj) appears. This term couples the density and orientation 
perturbations already in the linear order. As we will show in the following, this new 
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crosslink-induced anisotropic coupling modifies the density instability so it becomes 
transverse to the direction of polar orientation (in both the linear and nonlinear regime). 

Let us investigate the linear stability of the homogeneous polar solution of 
Eqs. (fT7|) and (fT8|) . describing a state with density po and polar orientation tq given 
by e(p — pc) = ^|Top. Without loss of generahty we set tq along x-direction, 
'''o = (l'''o|) 0). Linearizing the model equations around this state by making the ansatz 
{p, r^., Ty} = {po, Tq, 0} + {6p, 6tx, 6Ty} exp[a {k)t + ik^x + ikyu], one can deduce the linear 
growth rates cr as a function of the modulation wavenumbers kx,ky. For simplicity we 
set if = here. Finite but small values of H introduce a small drift but only slightly 
affect the growth rates. 

First consider the case without crosslinks {rj = 0). Then Eqs. f|T7|) . f|T8|) reduce 
to the model of Refs. [TTl [19]. There are three linear modes in the system. The 
two largest growth rates for long-wave perturbations are associated to a transverse 
orientational mode and to a mixed density-orientation mode. The third mode, related 
to the modulus of the orientation, is always damped. To leading order in kx, ky the 
transverse orientational mode reads 

a. = - (d., + (^n + + (20) 

and is thus always damped. For the mixed density mode one obtains 

= - [d, - ^) {kl + kl) . (21) 

Thus a density instability occurs at 

Po>P. = ^, (22) 

as already described in Refs. [ITJ [19], which to leading order is isotropic. Note that 
depending on the model parameters the density instability for rj = may also occur 
prior to the orientation instability, i.e. pd can be smaller than pc. 

A similar analysis can be done in the presence of crosslinks, rj = 1. While the 
orientational mode remains unchanged, for the mixed density mode one now obtains 

^^^^yi-y^^-^yi- (23) 

For perturbations in x-direction, i.e. parallel to the polar orientation, the density mode 
is damped. However, using the estimates from above, e ~ 0.273,^4 = 2.18, Dp = 1/32, 
and one obtains that the coefficient in front of ky is neg ative: Dp - < 0, 

i.e. transverse perturbations (i.e. with small k^ and finite ky) are unstable. 

Although this linear analysis reveals the possibility of a transverse instability in 
the presence of crosslinks, it is not clear if the density modulations perpendicular to 
the filament orientation really lead to bundle-like structures in the nonlinear regime. 
To investigate the long-term development of this instability, we performed numerical 
simulations of the full set of equations (fT7|) . (fTSjl . as described below. 
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Figure 3. Composite image of the density (color code: black low density, bright 
yellow high density) and the filament orientation field (arrows), a) The model of 
Ref. ^9J (?/ — 0) for po = 5 in the region of the isotropic density instability. Here 
the filament orientation is uncorrelated with the density gradient, b) The model with 
crosslinks (77 = 1) for a density of po = 6 displays pronounced bundles, and the local 
filament orientation is predominantly along the bundles. Other parameters values: 
if = 0.005,^2 = 0.6, C = 0.04. 



5. Numerical studies 

In order to study the system beyond the hnear regime, we performed numerical 
investigations of Eqs. ( |T7I) .( |T8l) . The studies were conducted in a SOL x SOL periodic 
domain, for different values of the parameter rj characterizing the concentration of 
crosslinks. Small amplitude noise was used as an initial condition for the r field, and 
P = Po + noise for the density field. Representative results for rj = 0,1 are presented 
in Fig. [3l In both situations, the simulations were performed in the regime where the 
homogeneous oriented state is unstable with respect to density fluctuations. However, 
depending on the value of the parameter rj, the manifestation of the instability is 
different. For r] = (without crosslinks), the numerical solution shows that the filament 
orientation and density gradients are mostly uncorrelated, cf. Fig. [3^. 

In contrast, for rj = 1 (with crosslinks), we observed that the instability 
indeed resulted in the formation of anisotropic bundles with the filaments' orientation 
predominantly along the bundles, as shown in Fig. [3]d. The bundles show a tendency 
to coarsen with time: small bundles coalesce into bigger bundles. The overall pattern 
is reminiscent of experimental observations of self-organization in both microtubules 
interacting with a mixture of motors of two different directions (kinesin and NCD) [12] 
and experiments on actomyosin where ATP-depleted myosin motors become crosslinks, 
cf. Fig.H 

In order to characterize the degree of alignment quantitatively, we calculated the 
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Figure 4. a) Composite image of the density (color code: black low density, red high 
density) and orientation field (arrows) for the model with crosslinks, same parameters 
as in Fig. [8)3 but at an early stage of evolution, b) Structure observed in microtubule- 
kinesin-NCD mixtures from Surrey et al. |12j . Here the two oppositely directed motors 
can be effectively mapped to the case of motor and crosslink, c) Experiment on 
actomyosin by Smith et al, where ATP-depleted oligomcric myosin- motors become 
crosslinks [20] . 



alignment coefficient between the density gradient Vp and the orientation r: 



where 0p and 0^ are the angles between Vp and the a;-axis and r and the x-axis 
correspondingly. The alignment coefficient C = 1 if the vectors Vp and r are everywhere 
perpendicular, and C = — 1 if they are parallel or antiparallel. 

We find an alignment coefficient of C = —0.0045 for the image shown in Fig. [3^ 
and corresponding to rj = (no crosslinks), con&ming that the fields r and Vp are 
practically uncorrelated. For the situation shown in Fig. [8)3 and corresponding to 77 = 1 
(crosslinks), we obtained a much larger value of C ~ 0.188, implying that the density 
gradient and the orientation are predominantly orthogonal. That means that density 
modulations are transverse to the orientation within a bundle. 

6. Conclusions 

As we have demonstrated above, the effect of crosslinks on the organization of polar 
filaments is twofold. First, the crosslinks, acting as hinges, allow zipping and result in 
the alignment of polar filaments by directional motion of molecular motors. Second, the 
ensuing polar state is unstable with respect to transverse density perturbations yielding 
bundles of oriented filaments, in contrast to the case without crosslinks in which the 
density instability is isotropic. 

This result has a simple physical interpretation. In the absence of crosslinks the 
motors tend to bring together the mid-point positions of microtubules, triggering an 
isotropic density instability. This instability is a direct counterpart of the aggregation 
or clustering in a gas of inelastic or sticky particles [23]. With a crosslink holding 
two ffiaments together at the intersection point, however, the motion of the ffiaments 
along the bisector is suppressed whereas the angular aggregation proceeds unopposed 




(24) 



Effects of crosslinks on motor-mediated filament organization 



12 



0.8 - 




0.6 - 

w 

0.4 - 
0.2 - 



0.2 0.4 0.6 0.8 1 

A(|)/7l 

Figure 5. The effective inelasticity factor e as a function of the initial angle between 
two rigid filaments A^/tt, see for detail Ref. [TOj . 

(furthermore, in fact it becomes much more effective). Thus crosshnks turn the isotropic 
instabihty into a transversal one. 

There are two experiments related to the model described here. The experiments 
on microtubules in the presence of two oppositely directed motor species, as reported 
in Ref. [I2], appear to produce the same qualitative result as the case of a single motor 
species mixed with crosslinks. This is because two motors moving from an initial 
intersection point in opposite directions along filaments also lead to their complete 
alignment. Additionally, our analysis possibly sheds new light on the interpretation of 
recent experiments by Smith et al. [20j on actin - myosin mixtures. In this experiment, 
no patterns were observed in a situation with abundant ATP. However, long dense 
bundles of actin filaments were observed when ATP was depleted by the multi-headed 
myosin motor constructs, for which it is known that in the absence of ATP they rigidly 
attach to the actin filaments and effectively become static crosslinks. Also in accordance 
with this interpretation, reinjection of ATP into the motor-filament solution resulted in 
a consequent dissolution of the bundles and homogenized the system anew. 

This experimental result can be interpreted as follows. As mentioned earlier, in the 
absence of crosslinks, the interaction of one motor with a filament pair does not result 
in complete alignment. In fact, the average decrease of the relative angle is of the order 
of 25-30% only, corresponding to a value of the restitution coefficient of 7 ~ 0.85 or 
to a value of the effective inelasticity factor e = 27 — 1 ^ 0.7. Since the inelasticity 
factor approaches 1 at large ■?/' (see Fig. |5]), it effectively produces a cutoff interaction 
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angle of the order of 0o ~ O.Gvr. Filament flexibility only slightly decreases this value 
[22]. Using the above values for 7 and (pQ, one finds from Eq. f|T9l) that the critical 
density needed for the orientational instability is about pco ~ 15.7. However, in the 
presence of crosslinks, the interaction becomes fully inelastic, and is described by the 
restitution coefficient 7 = 1/2. Also, the interaction leads to a complete alignment 
for any initial angle, so we can take (pQ = vr. The critical density for these conditions 
(7 = 1/2,00 = ^) is Pel ~ 3.66, which is more then four times smaller. Thus in 
the experiments, even if without crosslinks the motor density was not high enough to 
trigger the orientation transition, due to the crosslinking by ATP-depleted motors the 
system is likely driven beyond the threshold of orientation transition. Moreover, the 
oriented state is typically unstable with respect to a transverse instability leading to 
bundle formation, implying that bundles are competing with aster-like structures as the 
experimental pictures suggest. 

The inclusion of crosslinks in the model of filament interaction via molecular 
motors, Ref. [17], was straightforward and yielded nontrivial results. However, further 
generalizations of the model are needed. First, instead of the parameter rj interpolating 
between the cases with and without crosslinks, an additional field for the density of 
crosslinks should be introduced. In case of the actomyosin system, where ATP-depleted 
motors are acting like crosslinks, this field might be coupled via some simple reaction 
kinetics to the active motor density. Second, the role of filament fiexibility is worth 
investigating in some detail (cf. [22]). Furthermore, it is well known that in vivo, 
the cytoskeletal filaments are often met in a state of constant polymerization and 
depolymerization by means of ATP and GTP hydrolysis, another nonequilibrium process 
that is known to lead to structure formation [24l [25l [26] . The competition of the two 
main nonequilibrium processes in the cytoskeleton, active transport of the filaments 
by molecular motors and active polymerization of the filaments themselves, might lead 
to new and surprising behavior. Finally, an analysis of the homogeneous polar state 
in a filament-motor model with motor-induced drift, which we have neglected here, is 
adressed in Ref. [27]. 

We thank David Smith and Joseph Kas for stimulating discussions and for providing 
panel c) of Fig. HI This work was supported by the US DOE, grant DE-AC02- 
06CH11357. 



7. Appendix: Evaluation of the collision integral 

The first term of the collision integral, Eq. ([6]), can be simplified by integrating out the 

(5-function after having expressed 0i by (pi = 2(j) — 02 and ri by 

cos(?/'/2) , , ri J, , 

ri = rH ^ ' A r- ra) ' , A{r - ra) , (25) 

+ cos (■0/2) ri + cos^ip / 2) 

where = 0i — 02 as defined in the main text and where we have introduced the matrix 

^ , cos(20) sin(20) 

' sin(20) -cos(20) ' ■ ^ ' 
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(Note that after integrating over 5(0 — 0) the angle in the matrix A becomes 0.) Then 
one substitutes w = 2(0 — 02) and $ = Vi — r2. 

In the second term the 5-function leads to = 0i. After the suitable substitution 
w = — 02 this implies = — w/2. Finally one obtains the following simple form 



(27) 



X = I d^ I dwW{^, w) [P(r + Ai^, + w/2)P{r - ^2^, - w/2) 

-P{r,(f))P{r-^,(j)-w) 



with 



- 2cos(w/2)i-r7(i + /l) ^ 2cos(w/2)i + r/(i + A) 

^1 — 'a 7 77r\ , ^2 — 



(28) 



4cos(w;/2) ' 4cos(ti;/2) 

In case of = 0, i.e. in the absence of crosslinks, one regains = A2 = 1/2 as in the 
model of Ref. [IH]. 

To evaluate the spatial integral, one has to transform to the coordinates R = 
{R\\, R±) introduced in the kernel, Eq. fllip . These are connected to ^ via a simple 
rotation, 



R|l 
Ri 



Rj 



Rj 



cos (p sm I 
— sin cos I 



(29) 



and the collision integral becomes 



X 



dR^^dRj_ 



4>o 



dwW{Yi, w) P{r + AiR^i, + w/2)P{r - A^R^i, - w/2) 



P{r,ct))P{r-R^_^,2t^-w) 



(30) 
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